Stochastic Hard-Sphere Dynamics for Hydrodynamics of Non-Ideal Fluids

A novel stochastic fluid model is proposed with non-ideal structure factor consistent with compressibility, and adjustable transport coefficients. This Stochastic Hard Sphere Dynamics (SHSD) algorithm is a modification of the Direct Simulation Monte Carlo (DSMC) algorithm and has several computational advantages over event-driven hard-sphere molecular dynamics. Surprisingly, SHSD results in an equation of state and pair correlation function identical to that of a deterministic Hamiltonian system of penetrable spheres interacting with linear core pair potentials. The fluctuating hydrodynamic behavior of the SHSD fluid is verified for the Brownian motion of a nano-particle suspended in a compressible solvent.Comment: This work performed under the auspices of the U.S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344 (LLNL-JRNL-401745). To appear in Phys. Rev. Lett. 200

A Thermodynamically-Consistent Non-Ideal Stochastic Hard-Sphere Fluid

A grid-free variant of the Direct Simulation Monte Carlo (DSMC) method is proposed, named the Isotropic DSMC (I-DSMC) method, that is suitable for simulating dense fluid flows at molecular scales. The I-DSMC algorithm eliminates all grid artifacts from the traditional DSMC algorithm; it is Galilean invariant and microscopically isotropic. The stochastic collision rules in I-DSMC are modified to yield a non-ideal structure factor that gives consistent compressibility, as first proposed in [Phys. Rev. Lett. 101:075902 (2008)]. The resulting Stochastic Hard Sphere Dynamics (SHSD) fluid is empirically shown to be thermodynamically identical to a deterministic Hamiltonian system of penetrable spheres interacting with a linear core pair potential, well-described by the hypernetted chain (HNC) approximation. We apply a stochastic Enskog kinetic theory for the SHSD fluid to obtain estimates for the transport coefficients that are in excellent agreement with particle simulations over a wide range of densities and collision rates. The fluctuating hydrodynamic behavior of the SHSD fluid is verified by comparing its dynamic structure factor against theory based on the Landau-Lifshitz Navier-Stokes equations. We also study the Brownian motion of a nano-particle suspended in an SHSD fluid and find a long-time power-law tail in its velocity autocorrelation function consistent with hydrodynamic theory and molecular dynamics calculations.Comment: 30 pages, revision adding some clarifications and a new figure. See also arXiv:0803.035

Stacking Entropy of Hard Sphere Crystals

Classical hard spheres crystallize at equilibrium at high enough density. Crystals made up of stackings of 2-dimensional hexagonal close-packed layers (e.g. fcc, hcp, etc.) differ in entropy by only about $10^{-3}k_B$ per sphere (all configurations are degenerate in energy). To readily resolve and study these small entropy differences, we have implemented two different multicanonical Monte Carlo algorithms that allow direct equilibration between crystals with different stacking sequences. Recent work had demonstrated that the fcc stacking has higher entropy than the hcp stacking. We have studied other stackings to demonstrate that the fcc stacking does indeed have the highest entropy of ALL possible stackings. The entropic interactions we could detect involve three, four and (although with less statistical certainty) five consecutive layers of spheres. These interlayer entropic interactions fall off in strength with increasing distance, as expected; this fall-off appears to be much slower near the melting density than at the maximum (close-packing) density. At maximum density the entropy difference between fcc and hcp stackings is $0.00115 +/- 0.00004 k_B$ per sphere, which is roughly 30% higher than the same quantity measured near the melting transition.Comment: 15 page

Stanford Aerospace Research Laboratory research overview

Over the last ten years, the Stanford Aerospace Robotics Laboratory (ARL) has developed a hardware facility in which a number of space robotics issues have been, and continue to be, addressed. This paper reviews two of the current ARL research areas: navigation and control of free flying space robots, and modelling and control of extremely flexible space structures. The ARL has designed and built several semi-autonomous free-flying robots that perform numerous tasks in a zero-gravity, drag-free, two-dimensional environment. It is envisioned that future generations of these robots will be part of a human-robot team, in which the robots will operate under the task-level commands of astronauts. To make this possible, the ARL has developed a graphical user interface (GUI) with an intuitive object-level motion-direction capability. Using this interface, the ARL has demonstrated autonomous navigation, intercept and capture of moving and spinning objects, object transport, multiple-robot cooperative manipulation, and simple assemblies from both free-flying and fixed bases. The ARL has also built a number of experimental test beds on which the modelling and control of flexible manipulators has been studied. Early ARL experiments in this arena demonstrated for the first time the capability to control the end-point position of both single-link and multi-link flexible manipulators using end-point sensing. Building on these accomplishments, the ARL has been able to control payloads with unknown dynamics at the end of a flexible manipulator, and to achieve high-performance control of a multi-link flexible manipulator

Comment on `Universal relation between the Kolmogorov-Sinai entropy and the thermodynamic entropy in simple liquids'

The intriguing relations between Kolmogorov-Sinai entropy and self diffusion coefficients and the excess (thermodynamic) entropy found by Dzugutov and collaborators do not appear to hold for hard sphere and hard disks systems.Comment: 1 page revte

The Semiclassical Coulomb Interaction

The semiclassical Coulomb excitation interaction is at times expressed in the Lorentz gauge in terms of the electromagnetic fields and a contribution from the scalar electric potential. We point out that the potential term can make spurious contributions to excitation cross sections, especially when the the decay of excited states is taken into account. We show that, through an appropriate gauge transformation, the excitation interaction can be expressed in terms of the electromagnetic fields alone.Comment: 12 pages. Phys. Rev. C, Rapid Communication, in pres

Liquid-Solid Phase Transition of the System with Two particles in a Rectangular Box

We study the statistical properties of two hard spheres in a two dimensional rectangular box. In this system, the relation like Van der Waals equation loop is obtained between the width of the box and the pressure working on side walls. The auto-correlation function of each particle's position is calculated numerically. By this calculation near the critical width, the time at which the correlation become zero gets longer according to the increase of the height of the box. Moreover, fast and slow relaxation processes like $\alpha$ and $\beta$ relaxations observed in supper cooled liquid are observed when the height of the box is sufficiently large. These relaxation processes are discussed with the probability distribution of relative position of two particles.Comment: 6 figure